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Abstract 

Recent results on effects of Bose-Einstein symmetrization in a sys- 
tem of independently produced particles are interpreted in terms of 
statistical physics. For a large class of distributions, the effective sizes 
of the system in momentum and in configuration space are shown to 
shrink when quantum interference is taken into account. 

1. In recent papers |[ | we worked out the implications of the assump- 
tion that the identical pions generated in multiple production processes have 
only the correlations due to Bose-Einstein statistics. This was taken to mean 
that observable distributions can be evaluated in two steps. First, for the 
unphysical case of distinguishable pions independence is assumed, i.e. the 
poissonian multiplicity distirbution 

v N 

p(°)(iV) = e-— (1) 
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and the density matrix for each multiplicity in the form of a product of single 
particle density matrices. In momentum representation, for N particles we 
have 

N 

p£W) = n^<?o, (2) 
i=i 

where qt denotes the momentum vector of particle i and q is the set of all N 
momenta^. At this stage the momentum distribution is given by the diagonal 
elements of the density matrix 

rt°\q)=p%\q,q) (3) 

and is normalized to unity 

fdqp${q) = l. (4) 

In the second step the distribution is symmetrized, i.e. the resulting momen- 
tum distribution is obtained according to the prescription 

<M?) = 4 E P (0) (<?p,<7p) = E^p (0 W), (5) 
iV . ppi p 

where the first sum runs over all permutations P, P' of the momenta. As 
seen from (^), the momentum distribution f2jv(<z) is not normalized to one 
any more. Consequently, the multiplicity distribution changes from ([]]) into 

P(N) = C^ J Q N (q)dq (6) 

where C is a normalization constant. 

We found for the generating function of the multiplicity distributinon [|TJ 
and for the generating functional for the multiparticle correlation functions 
in momentum space explicit expressions in terms of the eigenfunctions and 
eigenvalues of the single particle density matrix p(°\ In the present paper 
we describe a reinterpret at ion of this model and show how the results can be 
simply obtained by using standard statistical physics. This explains directly 
the physical meaning of the results obtained in |2|. 

lr The special case, when is a gaussian, has been studied already by a number of 
authors @-§ 
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2. The density operator corresponding to the density matrix p^(qi,q'i) 
can be expressed in terms of its eigenvectors and eigenvalues 

p^ = J2\n>\ n <n\. (7) 

n 

Our first remark is that this operator corresponds to a single particle canon- 
ical distribution, if we make the identification! 

An = |e-^, (8) 

where (5 — -r^, e n is the energy ascribed to the state n, and Z is the nor- 
malizing factor which ensures that the density operator has the trace equal 
one 

E^ = L (9) 

n 

In terms of statistical physics, Z is the single particle partition function: 

Z = ^e- fe . (10) 

n 

The sum is over all states of the particle, thus a given energy e n can occur 
more than once. The hamiltonian with eigenstates |n > and eigenvalues e n 
is, of course, 

H = ^2\n > e n < n\. (11) 

n 

For instance, from results of [l[| and 0] we find that for the Gaussian single 
particle density matrix in one dimension 



" <0W) = TSa^"" " (12) 

where A 2 and R 2 are positive constants constrained by the condition RA > 
1/2 following from Heisenberg's uncertainty relation, the Hamiltonian is that 
of a harmonic oscillator with 

mu = h^ 2 (13) 



2 We consider only states with A„ > 0. The states with A„ = correspond to infinite 
energy and do not play any role in our argument. 
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and 

O DA 1 

Uu = -kT\n Z~ . (14) 
2RA + 1 v ; 

Inversely, for any Hamiltonian with a known set of eigenfunctions Tp n {q) and 
eigenvalues e n one can construct the corresponding density matrix (q, q') 
using the formulae (|?p and (|^). This allows to obtain a fairly large class of 
explicitly solvable models. 

Let us consider now a set of iV indistinguishable particles, where the 
Hamiltonian is the sum over all the particles of single particle Hamiltonians 
(|TT1). We assume that all these particles have been produced independently! 
and that the probability of producing N particles is v N . Let us consider first 
the subset of particles in state \n >. Since the number of particles in this 
subset is not fixed, we use the grand-canonical ensemble. The probability of 
finding exactly N particles in state |n > is 

P n (N) = —v N e-P N ^. (15) 

The parameter z/, known as fugacity, is related to the chemical potential p 
by the formula 

v = (16) 

and to v by v = vjZ. The normalizing factor 2 n , known as the grand 
partition function, ensures that the sum of all the probabilities equals one: 



,N 



*. = =j— K . (17) 

N=0 1 ^ 



This formula makes sense only if the geometrical series is convergent, i.e. for 
p < e n . The grand partition function for the whole system is a product of 
the grand partition functions for the independent subsystems, thus it is 



*=n r 



' (18) 



— ve P tn 



The grand partition function contains much information about the system. 
In particular it can replace the generating function for calculation of the 



3 These assumptions correspond to the mean field approximation in statistical physics. 
Thermodynamic equilibrium is not assumed. For a discussion assuming thermodynamic 
cqulibrium cf. S. 
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multiplicity distributions. For instance, the average number of particles in 
state n is 

d(lnZ n ) _ 1 _ v\ n 

dfi e/3(en-M)-l l-u\ n { f 

The first form is the well-known Bose-Einstein distribution, the second is 
the result obtained in [I]. Thus the symmetrization according to formula (|^) 
corresponds to the replacement of the Boltzman multiplicity distribution by 
the Bose-Einstein one. For the probability of no particle in state n (N n = 0) 
we get from ([15]) 

PJO) = — = i- pe~ pen = — — . (20) 
K J Z n N n + 1 K J 

Two limiting cases are of interest here. When all the occupation numbers 
Ni are very small P(0) = ll-fn(O) ~ e~S^ n = is very small for large 
N. When almost all the particles are in the ground state (BE condensation), 
P(0) ~ (1 + N)^ 1 and this probability becomes much larger. Thus the 
probability of producing an event with no 7r°'s (centauro event 0) is greatly 
enhanced by symmetrization, when a significant fraction of particles is in 
the condensate. This happens when [i — > €q from below, or in terms of the 
eigenvalues of the density matrix when u\q — > 1 from below. 

3. The momentum distribution and correlation functions can be de- 
scribed along the same lines. It was shown in p| that they can all be ex- 
pressed in terms of one function of two variables L(q, q') given by 

L(q,q') = j:MqW n (q)T^ r . (21) 
In particular, we have 

uj(q) = L(q,q) (22) 
for the inclusive single particle distribution, and 

K 2 ( qi ,q 2 ) = \L( qi ,q 2 )\ 2 (23) 

for the two-particle correlation function. 

Using and (|T6[), Eq.fl2"ip can rewritten as 

L(q,q') = E^(ghW) . (ew -, ) _ 1 (24) 
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Thus, it is just the element of the density matrix for the pure state n averaged 
over the Bose-Einstein distribution. Replacing the Bose-Einstein weights by 
the Boltzmann weights (^) one recovers the unsymmetrized single particle 
density matrix p^ (q, q'). At this point let us observe that the unsymmetrized 
density matrix p^ (q, q') corresponds better than the symmetrized one to 
the physical intuition since it contains the parameters which are easier to 
intepret. This rises an interesting question, to what extent the apparent 
parameters of the system, as determined from measured ui{q) and K2(q,q'), 
are modified with respect to the original "physical" parameters given by 

p (0 W)- 

This question was explicitly solved for the case of the Gaussian density 
matrix (corresponding to the hamiltonian of a harmonic oscillator): the re- 
sulting distributions are narrower (both in momentum and in coordinate 
space) than the original, unsymmetrized distributions []2|, [TJ. It is not obvi- 
ous, however, how general this result is. Below we discuss this problem for 
the case of hamiltonians of the form 

H=f- + V( X ). (25) 

with V(x) — > +00 when \x\ — > 00, to ensure the convergence of the trace of 
the density matrix. 

To this end let us first observe that the Bose-Einstein weights, when 
compared with those of Boltzmann, enhance the contribution from states 
with lower energy in the sum (0). Thus, if the wave functions ip n {x) and 
ip n (q) broaden with increasing energies e n (decreasing A n ), the symmetrized 
distributions are narrower than the unsymmetrized ones, i.e. we recover 
the qualitative result obtained for the Gaussian density matrix. For such 
a situation to occur, it is sufficient for example that the potential V(x) = 
X\x\ a where A and a are positive constants^. Unfortunately, the necessary 
condition(s) are not so easy to determine. It is also not immediately obvious, 
what this condition implies for the shape of the density matrix p^(q,q')- 

This argument also holds for the distribution in configuration space, as 
determined from the correlation function ^2(91, §2)- In this case we argue 
that the distribution in the variable q\—q2 broadens when the symmetrization 



4 The authors thank V.Zakharov for calling their attention to the virial theorem, which 
implies that. 
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effects are introduced. To this end let us observe that if all eigenvalues X n 
were the same, the sums in fl2TD and fl24"D would be proportional to 5(q — q'). 



Thus L(qi 1 q 2 ) would be infinitely narrow in qi — q 2 , corresponding to an 
infinite volume for particle production. Introducing the Boltzmann weigths 
(H) cuts the contribution from large energy levels and enhances the lower 
ones. The result is broadening of L. Since, as we argued above, the Bose- 
Einstein weigths work even stronger in the same direction, they will broaden 
the distribution even more. 

We thus conclude that the narrowing of the momentum and configuration 
space distributions as a consequence of the symmetrization of the multiparti- 
cle wave functions (first observed for the Gaussian density matrix) is actually 
a much more general phenomenon, valid for a rather broad class of distri- 
butions. More work is needed, however, to determine precisely the physical 
conditions which determine the behaviour of the effective widths of the dis- 
tributions after symmetrization. 

4. For fermions the difference is that P n (N) = for N > 1. Thus the 
formulae (|l7l) and ([Tj|) are replaced by 

Z& = 1 + j? e -/fc» ; = I](l + ve-^") (26) 

n 

The formula for the average population of state n becomes 

- kT " = -57 : = V 27 

dfx e^-ri + 1 1 + v\ n y ' 

and, of course, there is no condensation. After symmetrization the elements 
of the density matrix are replaced by 

LM = EMvm<l') * ) + 1 (28) 

n 

One sees from (E8|) that the Fermi-Dirac weights reduce the contribution 
from the low-energy states as compared to the Boltzmann weights which 
would be there without (anti) symmetrization of the wave function. Thus 
the argument given in the previous section implies that (anti)symmetrization 
of the wave function leads to broadening of the fermion distribution both in 
momentum and configuration space. 
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5. In conclusion, we have shown that the assumption that the identi- 
cal bosons (fermions) generated in multiple production processes have only 
the correlations due to Bose-Einstein (Fermi-Dirac) statistics can be refor- 
mulated in the language of statistical physics. It corresponds to the stan- 
dard Bose-Einstein (Fermi-Dirac) distribution of particles in the mean field 
approximation. Using this formulation we have shown for a large class of 
distributions that symmetrization of the multiparticle wave function implies 
narrowing of the particle spectra both in momentum and in configuration 
space (broadening is expected in case of fermions). This generalizes the re- 
sult found earlier for the gaussian distributions 
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